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Technical  unmary 




ABSTRACT 


Cii I 


A mathematical  model  for  nonlinear  heat  flow  in  a rigid  body  of 
material  with  memory  leads  to  the  integrodifferential  equation  problem: 


(HF) 


u (t,x)  = / a (t  - t)o(u  (t,x))  dt  + f (t,x) 
Z ■ 0 xx 


(0  < t < 00,  X c ]R) 


u(0,x)  = uQ(x) 


(X  e -R)  , 


which  is  analyzed  by  an  energy  method  developed  jointly  with  C.  M.  Dafermos 
Global  existence,  uniqueness,  boundedness  and  the  decay  of  smooth  solutions 
as  t -+•  “ are  established  for  sufficiently  smooth  and  "small"  data,  under 
physically  reasonable  assumptions. 
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Key  Words:  Nonlinear  Volterra  integrodiff arential  equations  (of  hyperbolic 
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Boundedness,  Asymptotic  behaviour.  Energy  method.  Resolvent 
kernels,  Frequency  domain  method.  Nonlinear  heat  flow  in 
materials  with  memory 


Work  Unit  Number  1 (Applied  Analysis) 


SIGNIFICANCE  AND  EXPLANATION 


Problem  (HF)  stated  in  the  Abstract  is  a mathematical  model  for  nonlinear 
heat  flow  (with  a finite  speed  of  propagation)  in  an  unbounded  one-dimensional 
rigid  body  of  material  with  "memory".  With  an  appropriate  interpretation  of 
the  term  o (u^)  x (HF)  has  a valid  physical  meaning  in  two  and  three  space 
dimensions . 

To  motivate  the  rather  technical  assumptions  concerning  the  kernel  a 
used  in  the  analysis  we  derive  (HF)  from  physical  principles  using  the  internal 
energy  and  the  heat  flux  in  such  a body,  expressed  as  certain  functionals  of 
the  temperature  and  of  the  gradient  of  the  temperature  respectively,  and  then 
applying  the  balance  of  heat  law.  The  assumptions  on  a are  then  motivated 
by  physical  and  partly  mathematical  arguments  involving  the  internal  energy 
and  heat  flux  response  functions.  Such  arguments  seem  to  have  been  omitted 
in  earlier  literature. 

Problem  (HF)  cannot  conveniently  be  solved  explicitly,  even  in  the  linear 
2 

case  (o(C)  = c £»  5 e 3R,  c a constant).  We  use  an  energy  method  developed 
jointly  with  C.  M.  Dafermos  (MRC  Technical  Summary  Report  #1876  - Communica- 
tions in  PDE,  to  appear)  to  explain  that,  under  physically  reasonable  assumptions 
on  the  functions  a,  a,  f,  u , problem  (HF)  has  a unique,  global  smooth  solu- 
2 0 

txon  u e C ([0,°°)  * ]R)  , provided  the  data  u^  and  f are  sufficiently 

smooth  and  "small".  Moreover;  this  solution  has  a finite  speed  of  propagation 
and  possesses  certain  boundedness  and  decay  properties  as  t -*■  ■».  The 
restriction  to  small  data  is  needed  to  preclude  the  development  of  "shock" 
solutions  due  to  possible  crossing  of  characteristics  of  a related  nonlinear 
wave  equation.  An  interpretation  of  the  result  is  that  the  memory  term  in 
(HF)  provides  a dissipative  mechanism  which  has  a smoothing  effect  for  "small", 
smooth  data,  under  physically  reasonable  conditions  on  the  kernel  a. 

The  energy  method  rests  on  an  auxiliary  inequality  (equation  (3.8)) 
which  is  derived  by  a simple  "frequency  domain"  argument  using  the  properties 
of  a.  The  remarks  following  Lemma  3.3  suggest  new  and  physically  meaningful 
conditions  on  a (expressible  in  terms  of  the  energy  and  heat  flux  response 
functions)  for  the  application  of  the  energy  method.  Such  generalizations 


summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A NONLINEAR  HYPERBOLIC  VOLTERRA  EQUATION  ARISING  IN  HEAT  FLOW 


John  A.  Nohel 


Introduction.  In  this  largely  expository  paper  which  is  based  on  recent  joint 

work  with  C.  M.  Dafermos  [1]  we  use  energy  methods  to  discuss  the  global  existence, 

2 

uniqueness , boundedness , and  decay  as  t «®  of  smooth  (C  ) solutions  of  the 
nonlinear  Cauchy  problem 

t 

ut(t,x)  = / a(t  - s)o(ux(s,x) J^ds  + f(t,x)  (0  < t < «,  x c H) 

<HF)  • 0 

u (0,x)  = uQ (x)  (x  e JR) 

for  appropriately  smooth  and  "small"  data  uQ,  f.  Here  a : [0,“>)  R,  a : R -+■  JR 
(o(0)  =0),  f : [0,«)  x r -*■  r,  uQ  : R -*■  R are  given  functions  satisfying  assump- 
tions motivated  partly  by  physical  considerations  sketched  below,  and  partly  by 
the  method  of  analysis;  subscripts  in  (HF)  denote  partial  derivatives.  Some 
comments  on  closely  related  initial-boundary  value  problems  are  made  following 
the  statement  of  the  main  result  (Theorem  2.1).  With  appropriate  interpretation 

of  the  term  <J(u  ) , problem  (HF)  has  a valid  physical  meaning  in  any  number  of 

2 

space  dimensions,  and  we  refer  to  [1,  Thm.  7.1]  to  such  a problem  in  R studied 
by  an  extension  of  this  method.  An  earlier  study  of  (HF)  by  Mac  Camy  [6]  is  based 
on  the  method  of  Riemann  invariants  and  is  therefore  restricted  to  a single  space 
dimension.  The  present  method  which  yields  more  widely  applicable  results  even 
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in  one  space  dimension  is  more  direct  and  simpler.  For  a similar  approach  to 
a problem  in  nonlinear  viscoelasticity  we  refer  the  reader  to  II,  Theorem  5.1] 
and  to  Mac  Camy  [71  for  the  Riemann  invariant  approach. 

To  motivate  the  assumptions  to  be  imposed,  particularly  with  regard  to 
the  kernel  a,  we  consider  briefly  the  problem  of  nonlinear  heat  flow  in  an 
unbounded  one-dimensional  rigid  body  of  a material  with  memory.  Let  u(t,x), 
e(t,x),  q (t , x) , and  h(t,x)  denote  respectively  the  temperature,  the  internal 
energy,  the  heat  flux,  and  the  external  heat  supply  at  time  t and  position  x. 
Following  Gurtin  and  Pipkin  [2] , and  also  Mac  Camy  [6]  , we  assume  a model  for  heat 
flow  in  which  e(t,x)  and  q(t,x)  are  respectively  the  following  functionals 
of  the  temperature  u and  of  the  gradient  of  temperature  ux: 

t 

(1.1)  e(t,x)  = bu(t,x)  + / B(t  - s)u(s,x)ds  (0  < t < =»,  x e R)  ; 

— 00 

t 

(1.2)  q (t , x)  = - / y(t  - s)a(u  (s,x))ds  (0  < t < ®,  x e JR)  ; 

— 00 

it  is  assumed  that  the  history  uQ(t,x)  of  u (and  hence  also  the  history 
u0x(t,x)  of  the  temperature  gradient)  is  prescribed  up  to  t = 0 and  for 
x e JR.  We  can  assume  without  loss  of  generality  that  uQ(t,x)  =0,  (t  < 0, 
x e JR) ; for  if  that  is  not  the  casi?,  it  is  easily  seen  from  what  follows  that 
this  merely  alters  the  forcing  term  f in  (HF) . It  is  reasonable  to  assume 
that  o i B + R is  a smooth  function  satisfying  o(0)  = 0,  o' (0)  > 0 (in  fact, 
o*(£)  e > 0,  (£  e JR)  - recall  that  for  linear  heat  flow  o(£)  = c£,  c > 0 
a constant) . We  shall  assume  that  b > 0 is  a given  constant  and  that  the 
given  smooth  "memory"  functions  &,Y  e L^(0,°°);  thus  e(t,x)  and  q(t,x)  are 
bounded  whenever  u(t,x)  and  ux(t,x)  are  bounded.  It  should  be  noted  that 
in  the  applied  literature  B,Y  are  linear  combinations  of  decaying  exponentials 
with  positive  coefficients. 

If  h (t,x)  denotes  the  external  heat  supply,  the  balance  of  heat  requires 

that 

(1.3)  e t (t , x)  = -q^tjx)  + h (t,x)  (0  < t < »,  x e jr)  . 

Substituting  (1.1),  (1.2)  into  (1.3),  and  using  the  assumption  that  u(t,x)  = 0 
for  t < 0,  x e JR  yields 
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g t t 

bu  (t,x)  + r—  / g(s)u(t  - s,x)ds  = / v(t  ~ s)o(u  (s,x))  ds  + h(t,x)  , 
* 9t  0 o xx 


or  equivalently 


t t 

(1.4)  bu  (t,x)  + / B(t  - s)u  (t  - s,x)ds  = / y(t  ~ s)a(u  (s,x))ds 

o z 0 X 


+ h(t,x)  - g(t)u(0,x)  , 


where  we  also  prescribe  the  value  u(0,x)  = uQ(x),  x £ 1R.  To  reduce  (1.4)  to 


(HF)  define  the  resolvent  kernel  p of  6 by  the  relation 


bp  (t)  + ( B* p ) (t)  = - 


(0  < t < «°)  , 


where  here  and  in  what  follows  * denotes  the  convolution  on  [0,t] . It  follows 


from  standard  theory  of  linear  Volterra  equations  that  (p)  has  a unique  solution 


P e C [0,°°).  If  g is  a given  function  on  (0,°°)  , the  solution  of  the 


Volterra  equation 


by(t)  + (B*y) (t)  = g(t) 


(0  < t < °°) 


is  given  by  the  variation  of  constants  formula 


y(t)  = 24—  + (p*9> 

b 


(0  < t < ®)  . 


Applying  (1.5)  with  y = ufc  and  g the  right-hand  side  of  (1.4),  one  sees 


that  (1.4)  is  equivalent  to  (HF)  with 


a (t)  = - y(t)  + (p*Y) (t) 

D 


0 < t < ® 


(1.7)  f(t,x)  = — (h(t,x)  - B(t)u  (x))  + p*(h(t,x)  - B(t)u.(x)) 

b 0 0 


(0  < t.  < 00 , X C ]R)  . 


To  motivate  the  assumptions  to  be  imposed  on  a,  we  note  from  (1.1)  and 


(1.5)  that 


u(t,x)  = - + (p*c)  (t,x) 

b 


(0  < t < °°,  x e ]R)  , 


whem  p is  the  resolvent  of  B.  If  p e L'(0,®)  then  f.  bounded  implies 


u bounded.  But  since  we  assumed  B e L (0,®),  the  Paley-Wiener  theorem 


s 
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applied  to  equation  (p)  yields  that  p e L (0,">)  if  and  only  if 

8(s)  # -b  for  Re  s > 0 

00 

where  8(s)  = / exp(-st) B(t)dt.  But  the  internal  energy  e is  positive  and 
so  from  (1.1)  ° 

00 

b + / 3(t)dt  = b + 6(0)  > 0 . 

0 

Therefore,  the  Paley-Wiener  condition  can  be  modified  to  the  statement: 
p e L1(0,“)  if  and  only  if 


Re  8(s)  + b > 0, 


(Re  s > 0) 


The  assumption  y e L1  (0 ,«)  *">  C*(0,<»),  together  with  (1.6),  then  imply  that 
a e (0,°°)  JO,00)  , and  an  easy  calculation  using  (p)  and  (1.6)  shows  that 


a (s)  = 


(Re  s > 0)  . 


b + 8(s) 


For  physical  reasons  and  (1.2)  one  needs  to  require  / y(t)dt  > 0,  and  so  it 
also  follows  from  (1.9)  that  ° 


(1.10) 


«,  / Y (t)  dt 

/ a(t)dt  = a (0)  = — - 

0 b + / 6(t)dt 


Physically  the  function  y represents  the  heat  flux  relaxation  function,  and 

it  is  reasonable  to  assume  that  y(0)  > 0.  It  then  follows  from  (1.6)  that 

a(0)  = y (0)  > 0.  If,  as  is  reasonable,  it  is  also  assumed  that  y'(0)  < 0 
b i 

and  8(0)  > 0,  one  also  has  from  (1.6)  that  a' (0)  = — y' (0)  + p(0)y(0)  = 

iy.(0)  _ SWiJCiO!  < 0. 

b b 

To  summarize,  the  following  assumptions  concerning  the  kernel  a in  (HF) 
are  reasonable  for  the  heat  flow  problem: 

oe 

(1.11)  a e C1  [0,«)  D ^(0,®) , a(0)  > 0,  a'  (0)  <0 , / a(t)dt  > 0 ? 

0 


as  we  shall  see  below  we  shall  require  additional  smoothness  of  a,  as  well  as 
a positivity  "frequency  domain"  condition  involving  the  Laplace  transform  of  a. 
The  implications  of  this  condition  are  discussed  in  Lemma  3.3  and  Remarks 
following  it. 

In  the  analysis  of  (HF)  which  follows  we  shall  impose  other  technical 
assumptions  (see  assumptions  (a),  (a),  (£) , (uQ)  below).  To  motivate  our  result 
for  (HF)  assume  for  the  moment  that  aft)  = a (0)  > 0 for  t ^ 0.  Then  (HF) 
is  formally  equivalent  to  the  Cauchy  problem 

(W)  u = a(0)o(u  ) + f , u (0,x)  = u (x)  , u (0,x)  = u (x)  = f(0,x)  . 

tt  XX  t U w i 

If  o is  "genuinely  nonlinear"  (o"(£)  * 0,  5 e 50  , Lax  13]  has  shown  that 
(W)  fails  to  have  global  smooth  solutions  in  time  (even  if  f = 0)  , no  matter 
how  smooth  ond  takes  the  initial  data  due  to  the  development  of  "shocks"  (the 
first  derivatives  of  u develop  singularities  in  finite  time  due  to  the 
crossing  of  characteristics).  Note  that  aft)  = a(0)  is  excluded  by  (1.11). 

Nishida  [8]  has  shown  that  for  the  nonlinear  wave  equation  with  frictional 
damping 

(W  ) u + ou  = a (0) a (u  ) , u(0,x)  = u (x) , u (0,x)  - u (x),  a(0)  > 0 , 

Cl  LL  t XX  U L I 

the  dissipation  term  au^,  a > 0,  precludes  the  development  of  shocks  if  the 
initial  data  are  sufficiently  smooth  and  "small".  The  proof  rests  on  the  concept 
of  Riemann  invariant  and  is  restricted  to  one  space  dimension.  For  a generaliza- 
tion of  Nishida' s method  to  the  forced  equation  (W  ) we  refer  to  Nohel  [9] . As 

a 

will  be  seen  in  (3.3)  below,  (HF) , under  physically  reasonable  assumptions,  is 
equivalent  to  a variant  of  (W^)  with  an  additional  memory  term  which  makes  our 
result  for  (HF)  (Theorem  2.1)  plausible.  The  Nishida  approach  applied  to  (HF) 
(necessarily  restricted  to  one  space  dimension)  was  studied  by  Mac  Camy  [6]. 

2.  Statement  of  Results.  We  make  the  following  assumptions.  Concerning  o let 
(o)  a e C3(K)  , a(0)  = 0,  o' (0)  > 0 , 

the  first  for  technical  reasons,  the  others  on  physical  grounds  (recall  that  in 
the  linear  version  of  (HF)  afu^)  = ux* • Concerning  the  kernel  a assume 
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(a) 


(i)  a e S(3) [0,»)  , 

(ii)  a (0)  > 0,  a'  (0)  < 0 , 

(iii)  tja(m)  e L3 (0,  ®)  (j,m  = 0,1. 2,3)  , 

oo 

(iv)  letting  a (in)  = / e lnta(t)dt  , 

0 

(2.1)  Re  a(in)  >0  (n  e JO  » 

where  [0,“)  is  the  set  of  functions  with  bounded  continuous  derivatives 

up  to  and  including  order  m.  From  (1.11)  above  the  conditions  a e C^fO,®), 

a e 1/(0,®),  a (0)  > 0,  a*  (0)  < 0,  a(0)  = / a(t)dt>0  are  reasonable  on  physical 

0 

grounds;  the  remaining  ones  are  needed  for  technical  reasons  of  the  analysis. 

See  additional  remarks  on  alternatives  to  the  frequency  domain  condition  (2.1) 
following  Lemma  3.3  below.  Concerning  the  forcing  term  f we  assume  (essentially 
for  technical  reasons) 


(f) 


f,f„,f  ,f  ,f  ,f  ,£.  e L2(0,“;L2(JO ) 

t X tt  tX  XX  ttt  ttx  txx 


The  initial  datum  uQ  is  assumed  to  satisfy 


V 


u„  , u^  ,u_  e L ( 3R)  . 

Ox  Oxx  Oxxx 


Note  that  in  (uQ)  no  assumption  is  explicitly  made  about  uQ(x);  however,  for 
the  particular  physical  problem  one  would  also  have  to  require  uQ  e L2(F)  in 
order  that  f defined  by  (1.7)  satisfy  (f) . Our  result  concerning  (HF)  is 
(see  [1;  Theorem  4.1]): 

2 

Theorem  2.1.  Let  the  assumptions  (o) , (a) , (f) , (uQ)  hold.  If  the  H (10 

norm  of  uQx  and  the  L2 ( [0,®) ;L2 (10 ) norms  of  f and  its  derivatives  listed 

in  (f)  are  sufficiently  small,  then  the  Cauchy  problem  (HF)  has  a unique  solu- 
2 

tion  u « C (10,®)  * R)  with  the  following  properties: 


(i) 

U. 

,u  ,u 

^,u^  ,u  ,ut. . ,u 

*u^  ,u  e 

L°° 

([0,«.);L2(»)) 

» 

t 

X 

tt  tx  XX  ttt 

ttx  txx  XXX 

(ii) 

u. 

,u  . , 

u ,u  ,u  ..  ,U.J. 

T 2 

,uJ.  ,u  e L 

([0 

,®);L2(K))  , 

t 

tt 

tx  XX  ttt  ttx 

txx  XXX 

(iii) 

ut 

(t,-) 

in 

L2(»)  as 

t ® , 

(iv) 

ut 

(t,x) 

,ux(t,x) ,utfc(t,x) 

•“tx11’10  '"x,"1 

,x) 

-*•  0 as  t -*■  ® 

uniformly  on 
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We  remark  that  conclusions  (iii),  (iv)  are  easy  consequences  of  (i)  , (ii) . 

It  also  follows  from  the  proof  that  the  solution  u has  a finite  speed  of 
propagation. 

We  also  note  that  the  results  of  Theorem  2.1  hold  (with  essentially  the 
same  proofs)  for  the  following  two  problems  of  heat  flow  in  a body  on  the 
interval  [0,1]  (see  [1,  Theorem  6.1]): 

(i)  (HF)  on  (0,®)  x (0,1)  with  homogeneous  Neumann  boundary  conditions 

at  x = 0 and  x = 1,  and  with  uQ(x)  prescribed  on  [0,1]; 

(ii)  (HF)  on  (0,®)  x (0,1)  with  homogeneous  Dirichlet  boundary  conditions 

at  x = 0 and  x = 1,  and  with  uQ(x)  prescribed  on  [0,1]. 

In  both  problems  assumptions  (a)  and  (a)  are  unchanged  while  assumptions  (f) 
and  (uQ)  hold  in  L [(0,1)  : L (1R)  ] and  in  l/(0,l)  respectively.  For 
problem  (ii)  one  adds  the  assumption  f(t,0)  = f(t,l)  = 0. 

For  a version  of  (HF)  in  two  space  dimensions  and  with  a similar  but 
technically  more  involved  proof  we  refer  the  reader  to  [1,  Theorem  7.1]. 

We  note  also  that  (HF)  is  of  the  abstract  form 

t 

u' (t)  + / a(t  - t ) Au ( t ) d t = f(t)  (0  < t < ®) 

0 

u(0)  = uQ  , 

£ 

where  A is  the  nonlinear  operator  Au  = - — o(u  ) plus  appropriate  conditions 

o X X 

at  ±®  or  suitable  boundary  conditions  at  x = 0 and  1.  Such  abstract 
problems  have  been  recently  studied  by  Londen  [4] , [5]  for  a class  of  kernels 
a(-)  which  are  positive,  decreasing,  convex  on  [0,®)  and  which  satisfy  the 
condition  a' (0+)  = -®  which  is  crucial  for  his  technique.  In  addition,  the 
solution  obtained  by  Londen  is  not  sufficiently  regular,  and  no  comparable  decay 
results  are  obtained. 

Finally,  we  observe  that  a comparison  of  Theorem  2.1  and  of  its  proof  with 
the  results  and  method  of  proof  by  Mac  Camy  [6]  shows  that  our  approach  is  more 
direct,  not  restricted  to  one  space  dimension,  and  yields  more  general  results 
(see  additional  remarks  following  Lemma  3.3). 

3.  Outline  of  Proof  of  Theorem  2.1.  To  simplify  the  exposition  we  shall  assume 
that  f ~ 0 in  (HF) , and  we  refer  the  reader  to  [1]  for  the  technically  more 
involved  treatment  resulting  from  f * 0;  the  method  is  unaltered  by  this 
simplification. 
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u (t,x)  + a (0)  (k*u  ) (t,x)  = a (0) a (u  (t,x))  (0  < t < °°,  x e ]R) 

tt  OZ.  TL  XX 

(3.2)  - 

u(0,x)  = uQ(x),  ut(0,x)  = 0 (x  e J?)  . 

Another  important  equivalent  form  of  (HF)  with  f = 0 resulting  from  (3.2)  is 

u (t,x)  + a (0)k (0)u  (t,x)  » a (0) o (u  (t,x))  - a(0) (k' *u  ) (t,x) 

tt  t x.x  t 

(3.3)  • (0  < t < °°»  x e 3R)  , 

u(0,x)  = uQ(x)»  ut(0,x)  = 0 (x  e ]R)  . 

Since  a(0)k(0)  = -a' (0)  > 0,  (3.3)  suggests  the  dissipative  mechanism  induced 

by  the  memory  term  in  (HF)  and  the  relationship  with  the  damped  wave  equation 

(W  ) . We  remark  that  if  f * 0 in  (HF)  one  adds  the  forcing  term 
a 

*(t,x)  = ft(t,x)/a(0)  + (k*ffc) (t,x)  + a (0) k (t) f (0,x)  to  (3.2)  and  (3.3),  and  v 

one  replaces  the  zero  initial  condition  by  ut(0,x  = f(0,x). 
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The  proof  of  Theorem  2.1  is  carried  out  in  two  stages:  (i)  A suitable 
local  existence  and  uniqueness  result  is  established.  (ii)  A priori  estimates 
are  established  to  continue  the  local  solution;  these  will  at  the  same  time 
yield  conclusions  (i) , (ii)  of  Theorem  2.1. 

b.  Local  Theory.  We  shall  make  the  temporary  additional  assumption  concerning  o : 

* 

(a  ) there  exists  pQ  > 0 such  that  a' (£)  ^ pQ  > 0 (£  e K)  . 

* 

Proposition  3.1.  Let  the  assumptions  (o) , (o  ) , (uQ)  hold,  and  let 

k',  k"  e C [0 ,°°)  L1(0,“>).  Then  the  Cauchy  problem  (3.2)  (resp.  (3.3))  has  a 

2 

unique  solution  u e c ([0,T  ) * H)  on  a maximal  interval  J0,TQ)  * 1R, 

Tq  £ +«,  such  that  for  T e [0,TQ)  one  has 

• oo  2 

(i)  all  derivatives  of  u of  orders  one  to  three  inclusive  e L ( [0 , T] ; L (1R)); 

(ii)  i£  Tq  < °°,  then 

00 

/ [uf(t,x)  + u2(t,x)  + u2  (t,x)  + •••  + u2  (t,x)]dx  -*■  ® as  t -*•  T . 
t x tt  xxx  0 

• OO 


We  remark  that  the  property  of  finite  speed  of  propagation  of  solutions  of  (HF) 

is  an  easy  consequence  of  the  proof  of  Proposition  3.1. 

The  proof  uses  the  Banach  fixed  point  theorem.  Let  X(M,T)  be  the  set  of 
2 

functions  u e C ([0,T]  * 1R)  for  any  T > 0 such  that  u(0,x)  = uQ(x), 
ut(0,x)  = 0 and  such  that 

(i)  u ,u  ,u  ,...,u  e L ( [0,T] ; L2(R))  and 

t X Ct  xxx 

oo 

(ii)  sup  / [u2(t,x)  + u2(t,x)  + u2  (t,x)  + •••  + u2  (t,x)]dx  <_  M2  . 

(0,T)  -<=»  Z X XXX 

Note  that  X (M ,T)  is  not  empty  if  M is  sufficiently  large,  and  that  if 
u e X (M ,T) , then 

(iii)  sup  { |u  (t,x) | , |u  (t,x) | , ju  (t,x) | , |u  (t,x) | , |u  (t,x) | } _<  M . 

[O.TJxr  x tx  xx 

2 

Let  S be  the  map:  X(M,T)  -*■  C ([0,T]  x ]R)  which  carries  a function  v e X(M,T) 
into  the  solution  of  the  linear  Cauchy  problem  (see  (3.3)  for  motivation) 
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(3.4) 


t 

? 


u (t,x)  + k(0)ui.(t,x)  = a(0)(a'(v  (t,x))u  (t,x)  - (k'*v  )(t,x)] 

tt  t X XX  t 

(0  < t < T,  x e 1R) 

u(0,x)  = uQ(x),  ut(0,x)  = 0 . 


Clearly  a fixed  point  of  S will  be  a solution  of  (3.2)  (respectively  (3.3)). 

To  apply  the  Banach  fixed  point  theorem  to  the  map  S one  first  shows  (by  an 
energy  argument,  for  details  see  (1.  Lemma  3.1])  that  if  M is  sufficiently 
large  and  if  T is  sufficiently  small,  then  S .naps  X(M,T)  into  itself. 

One  next  equips  X(M,T)  with  the  metric 

1 

00  J 

p(u,u)  = max  {/  [(ut(t,x)  - ufc(t,x))2  + (ux(t,x)  - ux(t,x))  ]dx> 

[0  ,T]  -°° 

By  the  lower  semicontinuity  of  norms  under  weak  convergence  in  Banach  space, 

X (M,T)  becomes  a complete  metric  space.  One  then  shows  that  for  M suffi- 
ciently large  and  T sufficiently  small  the  map  S is  a strict  contraction  of 
X (M ,T)  and  the  proof  of  Proposition  3.1  is  completed  in  a standard  manner  (for 
details  see  [1,  Lemma  3.2]). 

If  o,  k,  uQ,  u^  are  smoother,  the  solution  becomes  smoother.  A precise 
regularity  result  which  is  needed  for  the  a priori  estimates  is 

Proposition  3.2.  Let  the  assumptions  of  Proposition  3.1  be  satisfied.  In  addi- 
tion, assume  that 

(3.5)  a c C4(®)  , un  € L2(]R)  . 

Oxxxx 

Then  the  solution  u of  Proposition  3.1  has  the  addition  property 


(3.6) 


for  every  T < TQ,  where  ro,TQ)  x is  the  maximal  interval  of  existence. 

For  the  proof  see  (1,  Theorem  3.2]. 

£.  A Priori  Estimates  and  Continuati  wi.  We  wish  to  show  that  the  maximal 

interval  [0,T  ) of  Proposition  3.1  is  in  fact  [0,“) . Recall  that  the  local 

0 * 

theory  assumes  that  (o  ) is  satisfied;  this  assumption  will  be  removed.  The  a 
priori  estimates  will  be  deduced  from  equation  (3.2)  above.  We  shall  restrict 
the  range  of  ux(t,x)  for  a local  solution  u to  the  set  on  which  o' (•)  > 0; 
choore  c0  > 0 such  that 


1 


10 


(3.7)  o’(w)  > P0  > 0,  w e [-c0.cQ]  . 

We  wish  to  show  that  there  exists  a constant  y > 0,  y < cQ,  depending  on 
PQ'  / |k ’ (t)  jdt,  max  { | a’ ( • ) | , | a" ( • ) | , | o’" ( • ) | } , but  not  on  T > 0 


such  that  if  the  local  solution  u of  (3.1)  satisfies 


* 

(y  ) 


sup  { |u  (t,x) | , |u  (t,x) | , |u  (t,x) I , |u  (t,x) I > _<  y , 
0<t<T,xe]R  x tx  xx 


then  certain  functionals  of  the  solution  u are  controllably  small  (i.e.  these 
functionals  can  be  made  arbitrarily  small  by  choosing  the  initial  data  suffi- 
ciently small  in  the  appropriate  H norms).  More  precisely,  the  result  of  the 

a priori  estimates  which  follow  is  that  if  the  assumptions  of  Theorem  2.1  hold 

2 

(with  f = 0)  and  if  the  H OR)  norm  of  u is  sufficiently  small,  then  for 

Ox  * 

as  long  as  the  local  solution  u of  (3.2)  satisfies  the  condition  (y  ) for 
y > 0 sufficiently  small,  the  condition 

**  r 2 2 2 2 

(y  ) J [u. (s,x)  + u (s,x)  + u.(s,x)  + — + u (s,x)]dx 
t x tt  xxx 

.OO 

s 00 

+ / / [u2(t,x)  + u2  (t,x)  + •••  + u2  (t,x)]dxdt  <_  y2  (0  s < T) 

0 -OO  ^ ^t  XXX 

**  It 

is  satisfied.  The  inequality  (y  ) in  turn  implies  that  condition  (y  ) holds 
and  the  cycle  closes  in  a standard  manner  using  Proposition  3.1.  Thus  the 
maximal  interval  of  existence  of  the  solution  u(t,x)  is  (0,°°)  * » and  (y  ) 

holds  for  0 _<  s < oo.  This  proves  properties  (i)  and  (ii)  and  Theorem  2.1. 

** 

To  establish  the  a priori  inequality  (y  ) we  shall  work  with  the  equivalent 
form  (3.2)  of  (HF)  and  we  require  some  additional  properties  of  the  resolvent 


kernel 

k. 

Lemma 

3.3.  Let  assumptions 

(a)  be  satisfied  and  let 

k be  the  resolvent  kernel 

of  a ’ 

defined  by  equation 

(k) . Then 

(i) 

k e B2[0,<*)  , k(0)  = - 

a' (0)  _ 

2 0 ' 

[a(0)P 

'ii) 

k (t)  = k + K (t)  , k 

00  OO 

= — — > 0;  K(m)  e L1  (0 ,*)  , m = 0,1,2  , 
a (0) 

(iii) 

for  every  T > 0 and 

for  every  v e l2[0,T] 

there  exists  a number 

y > 0 

such  that 
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T T 

(3.8)  / v(t)  (k*v) (t)dt  ^ Y / v2(t)dt  . 

0 G 0 

Remarks  on  Lemma  3.3.  We  refer  to  [6,  Lemma  3.1]  for  the  proof  of  Lemma  3.3 
and  to  (1,  Lemma  2.1]  for  some  comments  and  corrections  of  that  proof.  Here 
we  make  some  additional  comments  concerning  the  energy  inequality  (3.8)  which 
is  of  independent  interest.  If,  as  is  the  case  here  k'  e L (0,°°),  the 
inequality  (3.8)  is  derived  by  the  following  simple  argument  (see  the  method 
of  (10,  Theorem  1]).  Extend  k'  evenly  for  t < 0,  and  let 


vT(t)  = 


v (t)  if  t e [0,T] 


0 otherwise  . 


T T T 

/ v (t)  ~ (k*v) (t)dt  = k (0)  / v (t)dt  + / v(t) (k'*v) (t)dt 
0 « 0 0 


T . 1 1 

= k (0)  / V2(t)dt  + 7 / v (t)  / k*  (t  - T)v(T)dTdt 


OO  OD  w 

= k (0)  / V2(t)dt  + jf  vT(t)  / k'(t  - T)vT(T)drdt  . 


Letting  v (n)  = / e intvT(t)dt,  (n  e TR)  , the  Parseval  and  convolution 

— <x> 

theorems  give 

rp  oo  oo 

/ V(t)  ~ (k*v)  (t)dt  = / |v  (n)  |2dn  + I |vT(n)  I2k' (n)dn  . 

—00 


4tt  ' T 

— oo 


But  k'(n)  = 2Re  k’(ih),  where  * is  the  Laplace  transform,  and  Re  k' (in)  = 

/s 

Re  [ink  (in)  - k (0) ] . Therefore, 

•p  OO 

f v(t)  (k#v) (t)dt  * —■  / |vT(n) |2Re[ink(in)]dn  • 

0 

Now  an  easy  calculation  from  equation  (k)  yields 


Re  ink (in)  = Re 


1 _ Re  a (in) 

a (in)  |a(in) |2 
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Thus,  to  prove  the  ineqaulity  (3.8)  it  suffices  to  establish  the  existence  of 
Y > 0 such  that 


Re  a (in)  > > 

|a(in) |2 


0 


(n  e »)  • 


This  is  precisely  what  is  done  with  the  aid  of  assumptions  (a)  in  [6,  Lemma  3.1], 
although  Mac  Camy ' s derivation  of  (3.8)  is  different  from  the  above  and 
unnecessarily  complicated.  The  above  suggests  that  the  frequency  domain  condi- 
tion (2.1)  in  assumptions  (a)  should  be  replaced  by  the  condition 


(S)  there  exists  o > 0 such  that  a Re  a(in)  >_  | a ( i n ) | 2 (n  c R)  . 

The  importance  of  condition  (S)  was  first  recognized  by  0.  J.  Staffans  [11] 
in  a different  context.  He  showed  [11,  Theorem  2]  that  condition  (S)  is 
satisfied  for  at  least  two  classes  of  kernels  of  importance  for  the  problem  (HF) 
* 1 

a e L (O,00)  n BV[0,®)  and  a strongly  positive  on  [0,®) , 

(i)  ‘ 2 -1 
i.e.  there  exists  an  e > 0 such  that  Re  a (in)  > e(l  + n ) > 


and 


(ii) 


a e 

on 


L^(0,°°)  and 
(0,®)  (here 


a and  -a'  are  nonnegative  and  convex 
a(0+)  = -a' (0+)  = +®  are  allowed)  . 


Staffans  also  gives  an  example  of  a kernel  which  is  a positive  definite  measure 
u satisfying  (S) , but  such  that  u is  not  strictly  positive  (Re  y(in)  > 0). 

Incidentally,  it  is  not  hard  to  show  that  if  a e L^fO,00),  a(0)  > 0,  and 
a is  either  strongly  positive  on  [0,®)  or  a satisfies  condition  (S) , then 
a'(0)  < 0.  It  is  also  important  to  notice  that  if  a satisfies  (S)  and  a(in) 
is  defined  as  a function  (e.g.  if  a € L1(0,00)),  then  a(in)  can  vanish  at 
most  on  a set  of  measure  zero  on  the  imaginary  axis. 

The  above  considerations  suggest  that  the  energy  inequality  (3.8)  is  true 
under  other  useful  conditions  which  are  much  more  general  than  assumptions  (a) , 
and  such  results  are  now  being  obtained. 

The  remainder  of  this  section  is  devoted  to  the  derivation  of  the  a priori 

** 

estimates  which  imply  (u  ).  Define 

w p 

(3.9)  W(w)  « / 5(t)d(  >yw^  we  I-c0,c0]  , 
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- 


where  the  inequality  follows  from  (3.7) . Let  u be  a local  solution  of  (3.2) 

* 

satisfying  (u  ) for  some  T * 0 and  0 < y < cQ.  Multiply  (3.2)  by  ufc  and 
integrate  over  [0,s]  * ]R.  Using  (3.9)  and  Lemma  3.3  (iii)  one  obtains  the 
estimate  (recall  we  are  doing  the  special  case  f = 0 in  (HP)  so  that  ufc(0,x)  = 0) 

00  OO  S 0° 

(3.10)  ^ / u 2(s,x)dx  + a(0)  / W(u  (s,x))dx  + y J J u^(t,x)dxdt 

2 -OO  t -CO  X 0 -00 


< a(0)  / W(u  (0,x))dx  (0  < s < T)  , 

— J x — 

— 00 

where  an  integration  by  parts  with  respect  to  x was  carried  out  in  the  term 


S oo 


/ / o(ux<t,x) )xut(t,x)dxdt  , 


0 -« 

followed  by  an  application  of  (3.9).  It  follows  from  (3.9),  (3.10)  that 

00  OO  S oo 

(3.11)  / u2(s,x)dx  + a(0)p  / u2(s,x)dx  +2 y J / u2(t,x)dxdt 

-00  t 0 -oo  x 0 Z 


< 2a  (0)  / W(u  (x) ) dx 
— 1 Ox 


(0  < s < T)  . 


Thus,  for  as  long  as  (y  ) holds  with  y < cQ,  the  quantities 

oo  oo  s 00 

/ u2(s,x)dx,  / u2(s,x)dx,  / / u2(t,x)dxdt 

— oo  —oo  Q — oo 

are  controllably  small,  uniformly  on  [0,T) . 

We  now  derive  two  additional  estimates,  the  first  by  differentiating  (3.2) 
with  respect  to  t,  multiplying  by  utt(t,x)  and  integrating  over  [0,s]  * R, 
the  second  by  differentiating  (3.2)  with  respect  to  x,  multiplying  by  u^tt.x) 
and  integrating  over  [0,s]  x Following  the  procedure  in  obtaining  (3.10), 

and  noting  that  since  ut<0,x)  = 0 one  now  has 

a2  8 

— - (k*ut)  (t , x)  = (k*utt)  (t,x)  , 

3t^ 

we  obtain 
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00  OO 

(3.12)  ^ / u2  (s,x)dx  + / o'(u  (s,x))u2  (s,x)dx 


+ a (0) y / / u (t ,x) dxdt  < ^ / u^t<0,x)dx 

0 — OO  — OO 


+ ^-r — / / o"(u  (t,x))u^  (t,x)dxdt 

£,  r\  ^ tX 

0 —00 


(0  < s < T)  , 


(3.13)  \ J u2  (s,x)dx  + / o'(u  (s,x)u2  (s,x)dx 

2 tx  2 J x xx 


+ a (0) y / / u2^(t,x)dxdt  <_  — / a ' (uQx (x) ) (x)dx 


0 -°° 


S oo 

+ f J a"  (u  (t,x))u  (t,x)u2  (t,x)dxdt  (0  <_  s < T)  . 

^ _ x tx  XX 

0 

We  add  up  (3.12),  (3.13)  and  we  claim  that  in  the  resulting  inequality,  and  as 
* 

long  as  (u  ) holds  with  y sufficiently  small,  each  term  on  the  right-hand  side 
is  either  controllably  small  or  can  be  majorized  by  the  sum  of  such  a quantity 
and  a quantity  that  is  dominated  by  one  of  the  dissipation  terms.  Indeed, 


since  from  (3.3) 


(3.14) 


utt(0,x)  = a(0)a(uQx)x  , 


the  L (K)  norm  of  utt_(0,x)  is  controllably  small.  The  two  space-time 
integrals  in  (3.12),  (3.13)  are  majorized  as  follows 


(3.15) 


a ~ 

,—  / / a"  (ux  (t,x) ) u^_x  (t ,x)dxdt  | 

0 -00 


S «> 

<_  ij-  max  |o”M  | / / u.  (t,x)dxdt  , 

r-Vco]  0 - 


]S”  7 

(3.16)  \-  J j o" (ux (t ,x) )utx<t,x)uxx  (t ,x)dxdt 


S oo 

< | max  | o **  C • ) | / / uxx<t,x)dxdt  . 

\c0.cQ]  0 — 
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To  estimate  the  integral  on  the  right-hand  side  of  (3.16)  we  have  from  (3.3) 

(3.17)  a(0)a'(u  (t,x))u  (t,x)  = u (t,x)  + a (0)k (0)u. (t,x)  + a(0) (k'*u  ) (t,x)  , 

x xx  tt  t t 

and  this  yields  from  (3.7)  and  standard  estimates 
s oo  s 00 

(3.18)  a2(0)p2  / / u2^(t,x)dxdt  _<  4 / / u2t(t.x)dxdt 

0 — 00  0 — OO 


+ 4a2(0)k2(0)  / 
0 


oo  00  s 

/ U2 (t,x)dxdt + 4a2 (0) (/  |k'(t)|dt)2  / 
r 0 0 


OD 

/ u2(t,x)dxdt 

— OO 


The  restrictions  imposed  on  p are  expressed  in  terms  of  parameters  fixed 
a priori.  For  example,  the  estimate  (3.15)  and  the  desire  to  absorb  this  term 
into  the  corresponding  dissipation  term  in  (3.13),  impose  the  restriction 
y max  | a " ( - ) ( <_  . The  combination  of  (3.12)-(3.18)  then  yields  that 


the  quantities 


00 

/ utt(s,x)dx, 

— OO 


oo 

/ utx(s,x)dx, 
— 00 


oo  g ao 

/ u (s,x)dx,  / / u ( t , x ) dxd  t , 
_ xl(  n _ tt 


s - 2 s °°  2 

/ / u (t,x)dxdt,  and  / / u (t,x)dxdt  are  controllably  small  uniformly 

OtX  XX 

—oo  0 — oo 

* 

on  [0 ,T) , provided  (p  ) holds  for  p sufficiently  small. 

To  get  the  final  set  of  estimates  one  assumes  temporarily  that  the  addi- 
tional smoothness  assumption  (3.5)  of  Proposition  3.2  is  satisfied,  so  that  u 
satisfies  (3.6).  We  form  the  second  derivative  of  equation  (3.2)  first  with 
respect  t,  multiply  by  utfcfc(t,x)  and  then  integrate  over  [0,s]  * R, 
secondly  with  respect  to  t and  x,  multiply  by  uttx^t,X^  ' anc*  t*len 
integrate  over  [0,s]  * R.  Following  the  above  procedure  and  using  equation 
(3.3)  to  compute  u... (0,x)  and  u..  (0,x)  a tedious  but  straightforward 
calculation  (see  [1],  section  4)  shows  that  the  quantities  / uttt(s,x)dx, 

<**  00  00  g oo  _ OO  g 00 

/ uttx(s,x)dx'  J utxx(s,x)dx'  ^ uxxx(s,x)dx'  ^ utttdxdt'  ^ / uttxdxdfc* 

—00  —oo  —oo  0 — oo  0 —oo 

s 0°  2 s 00  2 

/ / u dxdt,  / / u dxdt  are  controllably  small  uniformly  on  (0,T) , 

OCXX  . XXX 

—00  0 — OO 

* 

provided  (p  ) holds  for  p sufficiently  small.  Moreover,  the  detailed  estimates 
show  they  depend  solely  on  parameters  which  do  not  involve  the  additional  assump- 
tion (3.5).  Therefore,  a simple  density  argument  can  be  used  to  remove  the 
extraneous  assumption  (3.5). 


-16- 


Combining  all  the  controllable  estimates  obtained  above  yields  the 
**  * 

inequality  (y  ),  for  any  local  solution  u satisfying  (y  ) for  y > 0 

sufficiently  small.  This  completes  the  sketch  of  the  proof  of  Theorem  2.1. 
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